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We explore quantitatively the conditions in which matter (protons, neutrons) formed in the Early 
Universe during a period which spanned 10-50 pLS. We obtain all chemical potentials implied by 
the present day baryon-to-photon ratio. The baryochemical potential reaches its minimal value 
at the beginning of hadronization of the quark Universe, where pLb = 1.1 db 0.25 eV. In the mixed 
hadron-quark phase, a significant hadron sector electric charge distillation is found. 
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In the standard big-bang model, the large primordial 
baryon and antibaryon abundance formed at hadroniza¬ 
tion of the deconfined quark-gluon plasma (QGP) disap¬ 
pears due to mutual annihilation, exposing a slight net 
baryon number observed today 0. Applying the knowl¬ 
edge of equations of state of hadronic matter derived from 
the study of high energy nuclear collisions we consider 
quantitatively this evolution epoch of the early Universe. 

Our objective is to quantify the values of the chemi¬ 
cal potentials required to generate the observed matter- 
antimatter asymmetry, and to use this to quantify the 
magnitude of electrical charge distillation occurring dur¬ 
ing hadronization. While we do not enter here into the 
exploration of the consequences of these two findings, we 
note that the values of the chemical potentials allow de¬ 
tailed study of hadron and lepton abundances during the 
epoch the universe cools from the hadronization tempera¬ 
ture to nucleosynthesis temperature. The dynamical pro¬ 
cess of distillation can at QGP hadronization enhance an 
initially small baryon-antibaryon asymmetry to the value 
we encounter in our region of the Universe. 

The annihilation period began after the phase trans¬ 
formation from the QGP to a hot hadronic gas (HG), 
roughly 20-30//s after the big bang (when the Uni¬ 
verse was at a temperature of ~ 170 MeV), and con¬ 
tinued until the density had diminished to the level that 
a “nucleon freeze-out” occurred at approximately 1 ms 
(T ~ 35 MeV). In this scenario, annihilation of antimat¬ 
ter was very complete, quantified by the result that the 
energy fraction in baryons and antibaryons in the Uni¬ 
verse dropped from ~ 10% when Universe was about 
40 pLB old to ~ lO”"^ when it was one second old. 

The observational evidence about the antimatter non¬ 
abundance in the Universe is supported by the highly 
homogeneous cosmic microwave background derived from 
the period of photon decoupling Q. This has been used 
to argue that the matter-antimatter domains on a scale 
smaller than the observable Universe are unlikely Q ; oth¬ 
ers see need for further experimental study to confirm this 
result 0. 

The current small value of the baryon-to-photon ra¬ 
tio is the result of this near complete annihilation of the 
large matter-antimatter abundance. Other than a (rela¬ 
tively) small increase in photons during nucleosynthesis 
and electron-ion recombination, the baryon and photon 


numbers should be preserved back to the period of anni¬ 
hilation. Gonsidering several observables, a range of riis 
established and we use here the latest WMAP result j^, 
r] = UBlriy = h.llo;^ ^ 10“^°. 

The importance of rj is that it allows to determine 
the value of entropy per baryon S/B in the Universe, 
which is conserved in adiabatic evolution. At present 
the entropy is dominated by photons, and nearly mass¬ 
less (decoupled) neutrinos. It is straightforward to com¬ 
pute the entropy densities of these species from the 
partition function, and then to convert rj to S/B us¬ 
ing the photon number density. We obtain a value of 
S/B = 8.0/ry = 1.3 ± 0.1 x 10^°, assuming a lower neu¬ 
trino than photon temperature (photons are reheated by 
e+e” annihilation), and counting only left/right-handed 
nent rinos / antineutrinos. 

Equations of State We argue that the Universe was 
in chemical and thermal equilibrium around the time 
of hadronization, i.e., when the quark-gluon Universe 
turned into a hadronic Universe. To compute the ther¬ 
modynamic properties of the QGP and HG phases, we 
study the partition functions In Zqqp and In Zhg as de¬ 
scribed in Ref. 0: we employ latest QGP equations of 
state 0, which model in detail with properties of quan¬ 
tum chromodynamics (QCD) at finite temperature and 
finite baryon density obtained in lattice QGD. This ap¬ 
proach involves quantum gases of quarks and effectively 
massive gluons with perturbative QGD corrections ap¬ 
plied, and a confining vacuum energy-pressure compo¬ 
nent B — 258 MeV fm“^. In the HG partition function, 
we sum partial gas contributions including all hadrons 
from Ref. 0 having mass less than 2 GeV, and apply 
finite volume corrections 0. 

Our use of partition functions assumes that local ther¬ 
modynamic equilibrium (LTE) exists. Gonsidering the 
particle spectra and yields measured at the Relativistic 
Heavy Ion Collider at Brookhaven National Laboratory 
(RHIC), it is observed that a thermalization timescale 
on the order of rth ~ 10“^^ s is present in the QGP 
at hadronization 0. The mechanisms for such a rapid 
thermalization are at present under investigation. We ex¬ 
pect this result to be valid qualitatively in the primordial 
QGP phase of matter. This then assures us of LTE be¬ 
ing present in the evolving Universe. The local chemical 
equilibrium (LCE) is also approached at RHIC, indicat- 
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ing that this condition also prevails in the early Universe. 

To apply these experimental results, we recall that the 
size of a flat fc = 0 Universe evolves as i? oc where 
R is the size scale factor of the Universe M- This scal¬ 
ing is valid if the energy and momentum are dominated 
by radiation. Furthermore, if the expansion is adiabatic 
and energy conserving, then R oc T~^. The thermal- 
ization timescale is roughly rth ~ 1/nav, with n the 
particle number density, a the cross section for (energy¬ 
exchanging) interactions, and v the mean velocity. Al¬ 
lowing for the change in relative velocity and a reduction 
in density, we can expect an increase in Tth as we cross 
from the relativistic QGP to the HG phase having strong 
non-relativistic components, such that Tth ~ 10“^^ s for 
the HG at T = 170 MeV. For a roughly constant value of 
av during the expansion of the HG, the thermalization 
timescale increases to Tth ~ 10“^^ s at T = 1 MeV. At 
this point, the Universe is already one second old, so our 
assumption of LTE (and also of LCE) has a large margin 
of error and is in our opinion fully justified throughout 
the period of interest. 

Chemical equilibration timescales are significantly 
longer, due to smaller cross sections, than thermalization 
timescales. When chemical equilibrium cannot be main¬ 
tained in an expanding Universe, we find particle yield 
freeze-out. Near the phase transformation from HG to 
QGP, chemical equilibrium for hadrons made of u, d, s 
quarks is established. Hadron abundance evolution and 
deviations from the local equilibrium have not yet been 
studied in great detail in the early Universe. Our esti¬ 
mates suggest that the abundance of hadrons can be seen 
as being fixed by LCE down to a few MeV. Note that 
annihilation depletion of baryon density ceases around 
temperature T ~ 35 MeV. 

In a system of non-interacting particles, the chemical 
potential of each species i is independent of the chemi¬ 
cal potentials of other species, resulting in a large number 
of free parameters. The many chemical particle interac¬ 
tions occurring in the QGP and HG phases, however, 
greatly reduce this number. 

First, in thermal equilibrium, photons assume the 
Planck distribution, implying a zero photon chemical po¬ 
tential; i.e., jj,^ = 0. Next, for any reaction viAi = 0, 
where Vi are the reaction equation coefficients of the 
chemical species Ai, chemical equilibrium occurs when 
Vi^i = 0, which follows from a minimization of the Gibbs 
free energy. Because reactions such as f + f ^ are 
allowed, where / and / are a fermion - antifermion pair, 
we immediately see that fj,f = —fJ-f whenever chemical 
and thermal equilibrium have been attained. 

Furthermore, when the system is chemically equili¬ 
brated with respect to weak interactions, we can write 
down the following relationships El: 

Me ~ d'l^e ~ Mm ~ Mi^^ “ Mr ~ Mi^t — ^Mij (1) 
M« = Md-A^/, = Md , (2) 


with the chemical equilibrium of hadrons being equal to 
the sum of the chemical potentials of their constituent 
quarks; e.g. for S° (uds): meo = Mu + Md + Ms = 3Md ~ 
A/r/, and the baryochemical potential is: 
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fib ^-^{nd +^J.u) = (3) 

Finally, if the experimentally-favored “large mixing an¬ 
gle” solution [2 is correcU the neutrino oscillations 
Ve ^ ^ Vt imply that [ 13 : = fiv, 

which reduces the number of independent chemical po¬ 
tentials to three. We choose these to be Md, fie, and 
We next seek to establish what values of our independent 
chemical potentials are required to generate the observed 
matter-antimatter asymmetry. 

Single Phase In a homogeneous (i.e., single phase) 
Universe, we seek to satisfy the following three criteria: 
i) Charge neutrality {Q = 0) is required to eliminate 
Goulomb energy. This implies that: 

niifii, T) = 0, (4) 

i 

where Qi and Ui are the charge and number density 
of species i, and the summation is carried out over all 
species present in the phase. 

a) Net lepton number equals net baryon number [L = B) 
is required in many baryogenesis scenarios. This implies 
that: 


UL - UB ='^{L^ - =0, (5) 

i 


where Li and Bi are the lepton and baryon numbers of 
species i. 

in) Constant entropy-per-baryon (S/B). This is the 
statement that the Universe evolves adiabatically, and 
is equivalent to: 


Y.i Bini{fii,T) 


( 6 ) 


where Si is the entropy density of species i. 

Eqs. 0)-® constitute a system of three coupled, non¬ 
linear equations of three unknowns {fid, fie, and for 
a given temperature. These equations were solved nu¬ 
merically using the Levenberg-Marquardt method 0 to 
obtain Fig.^ which shows the values of Md, Me, and Mi/ 
as function of the age of the Universe (the top axis shows 
the corresponding temperature). As the temperature de¬ 
creases, the value of Md approaches weighted one-third of 
the nucleon mass ((2 to„ — mp)/3 ~ 313.6 MeV, see hor¬ 
izontal line in Fig.®. This follows because the baryon 
(proton and neutron) partition functions are in the clas¬ 
sical Boltzmann regime at these temperatures. 

The error bars arise from WMAP uncertainty in the 
value of 77 . The chemical potentials required to gener¬ 
ate the current matter-antimatter asymmetry are smaller 
than 1 eV (horizontal line in Fig.®. Near the Universe 
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FIG. 3: The fraction of baryons in the HG and QGP during 
phase transformation. The parameter /hg is the fraction of 
total phase space occupied by the hadronic gas phase. 


FIG. 1: Chemical potentials Hd, /ie, and fj,,^ aronnd the time of 
the QGP-HG phase transformation. The error bars arise from 
the uncertainty in rj. Insert — expanded view around the 
phase transformation. Horizontal and vertical lines inserted 
to guide the eye. 
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FIG. 2: The hadronic energy content of the luminous Universe 
as a function of temperature assuming a constant entropy-per- 
baryon number of 1.3 x 10^°. 

hadronization temperature Th additional sensitivity in 
the value of chemical potentials arises due to the uncer¬ 
tainty in the value of hadronization temperature. With 
6% uncertainty in Th = 170 ± 10 MeV we obtain 15% 
uncertainty in the minimal value of the d-quark chemical 
potential. This error is dominating the uncertainty of the 
minimal value of the baryochemical potential, and thus 
/TbUin =~ 1.1 ±0.25 eV. 

Given the chemical potentials of all particles we can 
trace out all particle densities in the early Universe and 
can, for example, derive the hadronic energy content 



FIG. 4: Net charge (including leptons) per net baryon number 
in the HG and QGP during phase transformation. Horizontal 
line at zero inserted to guide the eye. 

in the luminous Universe as a function of temperature, 
shown in Figure [3 The fraction of energy in baryons 
and antibaryons is roughly 10% at the QGP-HG phase 
transformation, but rapidly vanishes, becoming signifi¬ 
cant again only when the Universe has cooled, and enters 
nucleon rest mass matter-dominated era. 

Phase Transformation: QGP to HG During the 
QGP to HG phase transformation, when both phases co¬ 
exist, the macroscopic conditions i. - iii. above are no 
longer valid individually within either the QGP or HG 
phases, and that the correct expressions must contain 
combinations of the two phases. We therefore parame¬ 
terize the partition function during the phase transforma¬ 
tion as In Ztot = /hg In ^hg ± (1 - /hg) In ^qgp, where 
the factor /hg represents the fraction of total phase space 
occupied by the HG phase. The correct expression anal¬ 
ogous to Eq. 0 is: 

(5 = 0 = Vqgp ± riQ^ Vrg 
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= 


(1 - /hg) + /hg '• 


(7) 


where the total volume Vtot is irrelevant to the solution. 
Analogous expressions can be derived for Eqs. and ©• 

These expressions were used to obtain Fig. |21 which 
shows the fraction of the total baryon number in the QGP 
and HG phases as a function of /hg- We note that the 
QGP phase has a higher baryon number density than the 
HG phase throughout the transformation. This value is 
nearly constant, with « 3. 

In Fig.n it was assumed that the value of /hg evolved 
linearly in time and that the duration of the phase trans¬ 
formation was Th = 10 (IS. In reality, these quantities 
are sensitive to the properties of the equations of state 
and the dynamics of the phase transformation. A more 
complete evaluation will require application of transport 
theory and is beyond the scope of the current work. Our 
value of r/i is an estimate discussed in |^. 

Figure^lshows the net charge per baryon in each phase 
as a function of /hg- At its onset the small region of 
HG phase takes on an initial positive charge density, 
which can be attributed to the proton-neutron bias to¬ 
ward positive charge. As a result, the QGP domain takes 
on a (initially tiny) negative charge density. Such dis¬ 
tilled dynamical asymmetry in particle yields was previ¬ 
ously explored for str ang eness separation and associated 
strangelet formation [iRIlq . 

Since the sign of the effect seen in Figure0]is the same 
across the entire hadronization region, the total charge of 
the remaining QGP domains is ever-increasingly negative 


and one would expect development of electromagnetic 
potential, which effectively alters the values of chemical 
potentials for charged species. It is evident that the pro¬ 
cess of charge distillation will have a feed-back effect on 
the QGP-HG transformation, and that flows of particles 
will occur that will alter the uniformly small net baryon 
density This can affect (during the phase transfor¬ 
mation) any local initial baryon-antibaryon asymmetry, 
and may also serve as a mechanism for generating mag¬ 
netic fields in the primordial Universe |l8l| . Quantitative 
evaluation of this baryon asymmetry enhancement effect 
entails methods of advanced transport theory beyond the 
scope of this work. 

We have determined the chemical potentials required 
to describe the matter-antimatter asymmetry in the 
Early Universe. We have shown that the non-zero chem¬ 
ical potentials result in a charge distillation during the 
phase transformation, with the QGP and HG receiving 
negative and positive charge densities, respectively. We 
note that a separation of baryons and antibaryons into 
domains could maintain a homogeneous zero charge den¬ 
sity Universe, a phenomenon which could, e.g., play a 
significant role in amplifying a pre-existent, much smaller 
net baryon yield, e.g. arising within the realm of the 
standard model. The evolution of the chemical poten¬ 
tials allows to trace out the particle abundance in the 
early universe prior to nucleosynthesis. 
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